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Abstract 

At the high densities present in the interior of neutron stars, the neutrons are condensed into 
the 3 i-2 superfluid phase. While this condensation has little impact on the equation of state, it can 
have an important role in determining the low-temperature energy-momentum transport properties. 
The spontaneous breaking of baryon number by the condensate gives rise to the familiar Goldstone 
boson, but in addition, the spontaneous breaking of rotational invariance by the condensate gives 
rise to three Goldstone bosons, in general, one for each broken generator of rotations. These 
Goldstone bosons, which couple to the Z°, provide a new mechanism for neutrino emission. Using 
a low-energy effective field theory to describe the dynamics of these Goldstone bosons we estimate 
the neutrino emissivity of dense neutron matter and show that their annihilation is the dominant 
energy-loss mechanism over a range of temperatures. 
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I. INTRODUCTION 



One of the greatest challenges facing nuclear physics today is to understand the behavior 
of nuclear matter away from nuclear matter density. This impacts not only our understand- 
ing of the interior of neutron stars, but also our search for the deconfined phase of QCD in 
relativistic heavy ion collisions. We have a good phenomenological description of the inter- 
action between nucleons that extends up to momenta much greater than Fermi momentum 
in nuclear matter density Q. However, when nuclear matter is significantly compressed 
this phenomenological description becomes unreliable as one moves toward the deconfined 
phase of QCD. In fact, a number of more "exotic"' phases have been proposed as the ground 
state of nuclear matter at high densities, such as the color-flavor locked (CFL) phase at 
asymptotically high densities. 

At relatively low densities we have a good description of nuclear interactions which are 
dominated by the attractive S-waves, with higher partial waves suppressed by powers of 
the typical momentum. As the density is increased the repulsive nature of the S-waves at 
higher momenta becomes important and at ~ 1.5 times nuclear matter density, p nm , the 
average 3 p2 interactions are the most attractive suggesting the formation of a 3 P2 neutron 
condensate [3]. This has been known for sometime, and considerable work has gone into 
determining the magnitude of the 3 P2 gap as a function of density with the most sophisticated 
nuclear potentials 

QUI, 

and also with effective low-energy potentials [5j. In this work we 
point out that since a 3 P2 condensate spontaneously breaks rotational invariance, there 
will be three Goldstone bosons (angulons). In addition, baryon number is also broken, as 
in any superfluid, leading to the existence of another, well known, Goldstone boson. These 
modes will dominate the low-energy, low-temperature properties of the system. In particular 
they provide an important new mechanism for neutrino emission that is not exponentially 
suppressed at temperatures below the critical temperature, T c . 

We can make a rough estimate of the size of the contribution coming from angulon 
annihilation into a neutrino pair using dimensional analysis 

£ = G 2 F T 9 , (1) 

where Gp is the Fermi constant and T the temperature. This estimate assumes that powers 
of dimensionful quantities like the Fermi momentum kp or the value of the gap A are not 
relevant. We will argue that this is indeed justified. The temperature dependence does not 
have the characteristic exponential suppression ~ e -2A/T f oun( j j n processes involving gapped 
fermions but it is one power of T higher than the electron-electron scattering contribution. 
On the other hand it is not suppressed by the low electron density present in /3-equilibrated 
matter. There are two caveats with the estimate in eq. (JTJ). First, the size of the coupling 
between angulons and the weak neutral gauge boson is not at all obvious, if in fact one exists. 
Second, the annihilation process is proportional to two powers of the angulon density n(T), 
and since the Bose distribution function depends on the energy E = vp of the angulon, the 
density scales as n(T) ~ T 3 /v 3 . If the angulon speed v is small, the number of angulons in 
a momentum interval is greatly enhanced, and therefore so is the emissivity. For this reason 
it is important to obtain an estimate of both v and the dependence of the emissivity on v. 
We will estimate these factors using an effective theory to organize our arguments. A true 
model independent calculation is, unfortunately, not possible at the moment. 
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A. The effective theory 



When neutrons condense in the 3 P 2 superfluid phase the order parameter is given by 

i(n T a 2 a j V k n) = A jk , (2) 

where n are neutron field operators, V J = V J — V J , and a 2 o-j acts in spin-space. As the 
neutrons are coupled together in the 3 P 2 state, the order parameter A jfc is a symmetric, 
traceless tensor. In general, a traceless, symmetric tensor is determined by two orthonormal 
frames that diagonalize its real and imaginary parts, and two of the eigenvalues (the third 
one follows from the tracelessness condition). Depending on the value of the eigenvalues 
different phases with different unbroken symmetries arise. It is the dynamics of the system 
that determines which one of those phases is the true ground state, symmetry arguments 
alone cannot determine it. Unfortunately, not much is known about the form of the order 
parameter in dense neutron matter. At temperatures close to the critical one, arguments 
based on the Landau- Ginsburg energy suggest that Aq is real, but the relative sizes of the 
eigenvalues are hard to predict pj. At zero temperature those arguments fail and much less 
is known. Due to this uncertainty we will consider here a particular choice of the order 
parameter that makes some calculations feasible. The order of magnitude of the emissivity 
likely will not depend on this choice. 

We consider the phase where the equilibrium value of the gap matrix Ao has eigenvalues 
equal to the cubic roots of the unity. That is, there is an orthonormal frame where it can 
be written as 



A = A . (3) 




Rotational invariance is completely spontaneously broken down to a discrete subgroup and 
as a consequence there are three exactly massless Goldstone bosons, one for each rotation 
generator. In addition, as A / spontaneously breaks baryon number, there is a fourth 
(and well-known) Goldstone boson. Nuclear forces do not conserve spin and orbital angular 
momentum separately, due to the tensor and spin-orbit forces. The smallness of the 3 P 2 
gap, as well as its small mixing with the 3 F 2 channel, suggests that the effective strength 
of those forces at the Fermi surface is small. For the sake of argument it will be convenient 
to consider the enlarged rotation symmetry group including independent spin and orbital 
rotations. This group is explicitly broken by the tensor and spin-orbit interactions down 
to the diagonal group of combined spin and orbital rotations. The formation of the gap 
further spontaneously breaks this group down to a discrete subgroup which depends on the 
particular form of the eigenvalues of Ao. 1 The analysis that follows does not depend on the 
particular discrete group left unbroken and we will not discuss it further. This situation is 
represented by the diagram: 



1 For special cases, as when two of the eigenvalues are degenerate, a 0(2) subgroup of rotations is left 
unbroken. 
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SU S (2) <g> S0 L (3) <g> 17(1) tens ° r/spin ° rbit : S0j(3) ® Z7(l) 



(nn)^O 



(nn)^O 



discrete > discrete 

tensor/spin— orbit 

where the horizontal arrows represent explicit breaking and the vertical arrows denote spon- 
taneous breaking due to the pairing. Disregarding for the moment the effect of the tensor 
and spin-orbit forces we expect seven Goldstone bosons. We write the order parameter as 
A = UA = U^sAq^l, where U = e l2< ^° is a phase and £s an d £l are orthogonal matrices. 
The transformation rules under phase (e l9 ), spin (Rs) or spatial (Rl) rotations are: 

U - e™ U , A - Rg{O s ) A R T L (6 L ) . (4) 

At low energies the system can be described by the most general Lagrangian containing the 
low energy degrees of freedom (the Goldstone bosons). The Lagrangian invariant under the 
full SU S (2) x SO L (3) x 17(1) group is 



1 ' Tr [<% A<9 A f ] - f 2 Tr[d i Ad i A t ] - w 1 d i A\ k d ] A kj 



8A 2 

C2 



+ iH v Z$(Udotf - d UU r ) + iHaZ^iIJ 1 (Ad Q A* - d Q AA^)\ + ■ • ■ , (5) 

where Z® and Z? are the time and spatial components of the Z° boson. Missing from eq. (0) 
are terms that break the non-diagonal part of the rotation group, terms that do not vanish 
for non-unitary A /A, as well as terms with more derivatives whose contribution to low 
energy observables are suppressed by powers of the typical energy divided by A. 

To the terms explicitly shown in eq. (jHJ) we have to add terms that break separate spin 
and orbital rotations. Since the spin-orbit interaction drives the formation of the gap, its 
strength is suppressed by ~ l/log(A//x). One of its effects is to give a mass to three out 
of the seven Goldstone bosons, specifically, to the ones corresponding to opposite spin and 
orbital rotations £s = £ L . The remaining four are strictly massless, as they correspond to 
the breaking of the exact rotation and baryon number symmetry. Even though the size of 
this mass term is suppressed, for small enough temperatures, the number of these pseudo- 
Goldstone bosons is exponentially suppressed, and we will discard them. For the range of 
temperatures where it is relevant (if any), this extra degrees of freedom could contribute to 
the Vv emissivity and our calculation should be considered as a lower bound. On the other 
hand, the effect of the symmetry breaking terms on the interactions should be suppressed 
and we disregard them in our order of magnitude estimate. 



B. Determination of Low-Energy Constants 

In order for the low-energy effective field theory discussed in the previous section to 
be predictive, the a "priori unknown coefficients that enter v, w, f, H v and Ha, must be 
determined from QCD. Such a matching is not possible at this point in time but we are 
helped by the fact that all that is needed is information about neutron interactions close to 
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the Fermi surface. Thus we can imagine integrating out modes away from the Fermi surface 
and obtaining an effective theory valid for small excitations around the Fermi sphere. It 
is known that, to leading order, the resulting theory is very simple and contains, besides a 
kinetic term with a modified Fermi speed, only the interactions leading to the formation of 
the gap 0] • That means that any underlying theory resulting in the same Fermi speed and 
gap will lead to the same low energy properties of the system. We choose then a particularly 
simple one L N = Lg + with 

L N S = nt (ido + Xj + f^n - g SP2 xf 3 (n T a 2 a k V l nyn T a 2 a^n , 

= C V Z° rfn + C A Z» n^a l n - g ZVv V^{\ - , (6) 

where the tensor Xij — \{$ik$jl + $u8jk ~ is a projector onto the 3 P 2 channel, \i 

is the neutron chemical potential, and gsp 2 is an effective strong coupling constant in the 
3 P 2 channel. Not much is known about the renormalization of the weak interactions by the 
process of integrating out modes away from the Fermi surface. In particular we do not know 
whether the form above is universal, in the same sense as the strong part of the Lagrangian 
is. We have kept the same form of the interaction as in the vacuum and also retained only 
the leading order terms in the derivative and multi-body expansion of the weak interactions. 
The multi-body operators are analogous to Lia m the pionless effective field theory (lp| . 

The fact that the strong interaction in eq. (jHJ) occurs in the 3 P 2 channel is suggested 
by the vacuum value of the phase shifts for nucleon-nucleon scattering, and is supported 
by sophisticated nuclear models. This is an assumption underlying our work. The neutron 
mass is renormalized when modes far from the Fermi surface are integrated out of the theory 
and thus M is the renormalized mass or the "in-medium" effective mass. An analogous 
renormalization occurs for the coupling to the weak currents and the chemical potential \i. 

The model in eq. © favors neutron spin-pairing and the formation of a gap in the 3 P 2 - 
channel. The strong coupling gzp 2 in eq. © is traded for the neutron gap, via (A )y = 
— gzp 2 (n T (—p) a 2 a.i pj n{p)). The tilde is to denote the model-dependence of this gap. The 
propagator for neutrons in this condensed phase is 

iSfr, P) = 2 - * 2 . ( , A P \ + €p , "^tfW \ f (7) 

where e p = |p| 2 /(2M) - fi. 

The neutral current couplings in eq. © are 

2 G F Ml 2 ~ 2 G F M\ 2 ~ 2 G F M 2 Z 

9zvv ~ ' v ~ v ^7T ' A " A ' [ } 

where Cy = —1, constrained by vector current conservation, and Ca = gA + As ~ 1.1 ± 
0.15 |3. g A ~ 1.26 is the nucleon isovector axial coupling that is well measured in nuclear 
/3-decay, while As is the matrix element of the strange axial-current in the proton that is 
measured in deep-inelastic scattering and neutrino-nucleon interactions. 



C. Estimate of fa and 

In order to determine the parameters for the phonon it is convenient to introduce 
a fictitious U(l) gauge-symmetry into the theory. We then require that the low-energy 
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FIG. 1: The one-loop diagrams that determine the decay constants /o and vq. The solid line 
denotes a neutron, while the wiggly line denotes the fictitious gauge field Aq. 

effective field theory reproduce matrix elements of the underlying theory. At leading order 
in the derivative expansion, this amounts to replacing partial derivatives with covariant 
derivatives: 

d^n -> D„n = d^n + lA^n , <9 M A -> D„A = d^N + HA^A , (9) 

where A^ is the fictitious gauge field associated with the fictitious U{1) gauge-symmetry. In 
the underlying theory A couples to the neutron density and acts as a chemical potential. 
The correlation function (AqAq) is the linear response function determining how the density 
changes due to a change in chemical potential. In other words, it is the density of states at the 
Fermi surface dN/dfj, = Mkp/n 2 (up to corrections of order ~ A///), where kp is the Fermi 
momentum. In a diagrammatic calculation, (AqAq) is given by the two diagrams of Fig. (fl|). 
Using the propagator in eq. (J7|) we find that the two diagrams give equal contributions 

, . . , , , d A k -(k + e k ) 2 + k.A Alk . f d 3 k k 2 A 2 



(2tt) 4 (k 2 -e 2 -k.A Al.k) 2 ~ J (27r)3( e 2 + A .2A2)l 
Mkp A 2 , 

*i—£ + Q(—), 10 

where we used the fact that the integral is dominated by momenta around the Fermi surface. 
Similarly, one can compute (AiAj). In the underlying theory, in Coulomb gauge V- A = 0, it 
is given by the same two diagrams in Fig. (JTJ plus a tadpole graph with a —AiA 1 rvn/{2M) 
vertex, as shown in Fig. (J2J). In contrast to the calculation of (A A ) there is a relative 
minus sign for the anomalous graph in Fig. (J2J) due to the derivative coupling at the vertices. 
Consequently, these two graphs cancel, leaving the tadpole contribution 

N k 3 

i(AiAj) = -i—Sn = -i — if— <5« . (11) 

In the effective theory (AqAq) and (AqAq) are given by the tree level contributions and 
we find 

Mk F 2 
i{A A ) = % — — = if , 

7T 

iiAAj) = -i^i^Sij = -iMSa . (12) 



/o 2 = ^ , v = ^. (13) 



We conclude that 

Mk F v\ 

Tl Z 3 

The breakdown of rotational invariance plays no role in the propagation of <ft and, conse- 
quently, the value of vo agrees with a general analysis which assumes rotational invariance 
(after a suitable generalization to the non-relativistic context) [l3| . 
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D. Estimate of / and v 

To estimate / and v we consider the fictitious approximate gauge spin rotation symmetry 
n — > e l 2" e n. The microscopic theory possesses this symmetry (up to corrections due to 
tensor /spin-orbit forces) if the neutrons couple to the fictitious gauge field through a 
covariant derivative 

D ll n = d ll n + iB)^n. (14) 
The effective theory should then be written in terms of covariant derivatives of the A field 

D^A — <9 M A + iB^A . (15) 

We then match the quantities (BqB 3 ) and (B^Bj) in the microscopic and the effective theory. 
The graphs contributing to them in the microscopic theory are the same ones as in the case 
of the A field correlators. We find 

■ / r* W \ 1 f d * k -(^o + e fc ) 2 Tr[ C rV]-(Ao.A;)a(AS.A;) 6 Tr[ ( TVVV'] 



4 7 (2tt)' {kl-elf 

• r Mk F 

where we used Tr[a k a % a l a^] = 2(5ikSji — SkiSij + SjkSu). Notice that the denominators of the 
fermion propagators are spherically symmetric, as k.\A \ 2 .k = k 2 A 2 . We also match the 
spatial part: 

1 , 2 f d 4 k -(k + e^Tia'a^-hk) + Ti(A .k) a a a a\Alk) b a b ^k k k l 



7R!R n _ 1 - 2 f d*k -(fco + e^Trtry^-fcfcfcQ + TrCAo 
H-Ofc-Dj/ - 4M2 ^ J (27f)4 {kl-ElY 



1 - r Tr[^]5 kl f d4k l(ko " ek) 



2M L '-J (2tt)' k 2 - e 2 k 

(17) 



4tt 2 M 



~^'*« + ^( A o,Ajl + Ao/Aji + A^Ajj + A^Ag) 



where the integrals were computed up to corrections of orders 0(A//x) but the matching 
is valid only up to much larger terms, of order g ~ l/log(A //i), since the spin rotation 
symmetry is only an approximate symmetry. 

The same matrix elements are given in the effective theory by 

4'2 I'2 2 -f'2 2 

i{BiBl) = iL.5* , i(BlBf) = -il^-^S^-i^iAliJ'P + Pj^Ao + c.c.UilS) 
and by matching to the expressions in the full theory using the relation 
(A^Ag + A^Aji + A^Ag + A^Ajj) = 4AW« - (Aj( JV' + J^)A + c.c.)w (19) 



gives 



/> = ^ . ^ = i4 , ^ = ±* ^ (20) 
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FIG. 2: The one-loop diagrams that determine / and the "speeds" v and w in eq. The solid 
line denotes a neutron, while the wiggly line denotes the fictitious gauge field B{. The solid circle 
denotes an insertion of the tadpole vertex. 

We parameterize the scalars fields asU = e 2 *^ , £ = e l ^ J ' 7T ^ so the fields and rc l will be 
canonically normalized. However, the space derivative terms mix the different components 
of 7r\ The matrix for the quadratic of the 7Tj Lagrangian reads 



G(p) 



1 /" d 4 p 



27 (2n)F i (~ p } G ^ ijiri ^' 



(21) 



4 PxPy 
w 2 

—PxPz 



pi - v 2 p 2 - \(p 2 x + pz< 

o 

w 

—PyPz 



\PxPz 
\VvPz 

pi - v 2 p 2 - \{p 2 x +p: 



where p 2 = p 2 +p 2 +p 2 . The kinetic part can be diagonalized by using the fields o 1 defined 
by 

7Ti(p) = K ia {p)a a . (22) 

The explicit form of K ia (p) can be easily found but we will omit it here since is not very 
enlightening. 



E. Estimate of the Weak Couplings Ha,H 



v 



Looking at the coupling of the fictitious fields and to the neutrons we see that we 
can identify 



A 
Bl 



-2C A (z°y 



from which we can determine Hv and Ha to be 



fo n 

T U V 



H t 



4A 2 



(23) 



(24) 



Among other contributions, these weak coefficients receive corrections of order ~ 1/ log(A//x) 
coming from the strong interactions that are not invariant under local spin rotations. In 
terms of the eigenmodes a 1 we have the couplings 

£zo = -f C v Z°d <p - fC A ^Z°K ia d a a + ^Z° k e^ k K m (p)K jb (k)a a (p)d a b (k) + ■■■ (25) 
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II. NEUTRINO EMISSIVITY 



The Goldstone bosons identified in this work will contribute to the neutrino emissivity of 
dense nuclear matter at finite temperature. The main mechanism is the annihilation process 
ajCtj — > Vv. The amplitude for the annihilation of a, and otj is given by 

M a6 = — y= [H abk E b (k) + H bak E a (p)} z/7 fc (l-7 5 )z/ ( no sum implied), (26) 

where Ei(p) is the energy of the angulon ai with three momentum p = \p\ and H abk = 
Y,ij=xZ ijk Ki a {p)K jb {k). ' 

The emissivity due to these annihilation processes is defined as the energy loss due to 
neutrino emissions per unit time per unit volume, and for a,<x,- — > Vv it is given by 

[ d 3 p u d z py (fp d 3 k fjpfhW 

tab J {2^f2E u {2-Kf2E v {2*f2E a {p) {2irf2E h {k) H a[P)) H b[ )] 

(E^E^) (27r)U^{p + k-p v -p T r) ^|M ab | 2 , (27) 

s,s' 

where the angulons a a , a b carry incoming-momenta p, k respectively, as indicated. n(E) = 
[exp(E/T) — is the Bose distribution function and ^ |M a f,| 2 is the spin summed squared 
matrix element obtained from eq. ()26|) . Using the Lorentz invariant quantity 

= / iSk pfe (2 * )4 * (4,( « - ^ - ^ 1 ~ 75)1 

= ^ (4.9, - <zV) %o)% 2 ) , (28) 

for the neutrino phase-space integration gives 

2 

' rk (E a (p) + E b (k)) n(E b (p)) n{E b {k)) 




(2ir) 3 2E a (p) (2n) 3 2E b (k) 

H abk E b {k) + H bak E a {p) } [ H abl E b {k) + H bal E a (p) } 
((p + k) fe (p + k), + (E a (p) + E b (k)) 2 8 kl -\p+ k\ 2 5 kl ) 

6({E a {p) + E b {k)) 2 -\p + k\ 2 ) . (29) 

The phase space integral is complicated by the angular dependence of the vertex tensor 
H a bk but we can extract its dependence on T and v (= w). First, redefine Ei(p) — vEi(p) 
in order to have the factors of v explicit. Then expand the step function as 6(v 2 (Ei(p) + 

^•(fc)) 2 -|p + A:| 2 )) ^ 6'(-|p + A:| 2 ))+^ 2 (^(p) + ^(fc)) 2 5( ? ; 2 (^(p) + ^(A:)) 2 -b + A:| 2 ))H . 

The first term does not contribute to the integral and we are left with 

£ „ v * f —„ n(vE a (p)) n(vE b (k)) ( ' E a (p) + E b {k) 

J (2^2E n (p) (2n) 3 2E b (k) V U JJ V K JJ \ W V ; 



(27r) 3 2£ a (p) (2^2E b (k) 
[ H abk E b {k) + H bak E a {p) ] [ H abl E b {k) + H bal E a (p) } 5(-\p + k\ 2 )) 

rp9 



V 3 



(30) 
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where in the last step we rescale the momenta as p — > Tx/v, k — > Ty/v and used the fact the 
delta function restricts the six-dimensional integral to back-to-back pairs with p = —k only. 
We have checked that in the analytically calculable case where the vertex is independent of 
the angle, the scaling above is obeyed. 

Our estimate for the annihilation process is then 



where we used for the numerical estimates v = 0.15 (corresponding to kp ~ 308 MeV, 
M = 940 MeV, neglecting the renormalization of the mass and chemical potential), T 9 = 
T/(10 9 K) and the factor of six is due to the six possible combination of angulon pairs 
annihilated. 

At typical temperatures, e.g. T ~ 3 x 10 8 -K" ~ T c /10, and densities the emissivity due to 
electron bremsstrahlung is of order £ e ~ 10 10 erg cm~ 3 s _s [14], which is significantly less 
than that due to angulon annihilation of 8, aa ~ 10 12 erg cm -3 s~ s , where we have used vf ~ 
0.33. Further, comparing with processes involving the neutrons near the Fermi surface, such 
as modified Urea and neutron bremsstrahlung, one finds that at temperatures much below 
the critical temperature, T c , where such processes are exponentially suppressed HQ, the 
annihilation of angulons is likely to dominate the emissivity as it is power-law suppressed 
only. 

III. CONCLUSIONS 

We have pointed out that in the 3 P2 neutron condensed phase that is favored for densities 
greater than ~ 1.5 p nm , there are Goldstone modes that contribute to the neutrino emissivity 
and the energy-momentum transport properties. Using effective field theory arguments we 
estimated the neutrino emissivity from finite temperature, superfluid neutron matter for a 
particularly simply form for the gap, one that gives a unitary order parameter. We showed 
that the emissivity can be related to the Fermi speed and the parameters determining 
the anisotropy of the gap, but is not strongly dependent on the value of the gap itself. 
By using reasonable estimates of these parameters we estimate an emissivity larger than 
other processes involving neutrons, which are exponentially suppressed, and larger than that 
from electron bremsstrahlung, for the densities and temperatures relevant to neutron stars. 
Further, these Goldstone bosons will likely dominate other low temperature observables such 
as neutrino opacity and viscosities. Our calculation was dependent on a particular choice of 
the form of the gap parameter. This highlights the fact that a determination of the actual 
phase of cold neutron matter would have a larger impact on understanding the cooling 
processes of the 3 P2 phase than a precise determination of the size of the gap. A better 
assessment of the impact of angulon annihilation in the cooling of neutron stars requires the 
rates computed in this paper be inserted into a realistic cooling code. 
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